The bound state problem of S-wave heavy quark meson-aitimeson systems 
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We investigated systematically whether the S-wave (Qq) meson and the (Qq) meson may form S- 
wave bound states in a chiral SU(3) quark model by solving the resonating group method equation. 
Here Q = c or 6 and q — u, d oi s. Our preliminary calculation disfavors the existence of / = i 
{Ql)-{Qs) molecules {I = u,d) while favors the existence of isoscalar BB, B* B* (J=2) and BB* 
(C=-|-) molecules. The existence of isovector (charm-anticharm) and (charm-bottom) molecules is 
also disfavored. Therefore the resonance-like structure Z^(4051) is unlikely to be an S-wave D* D* 
molecule. 

PACS numbers: 12.39.-x, 12.40.Yx, 13.75.Lb 



I. INTRODUCTION 

The molecular picture was widely used in discussing 
the strange states, such as /o(980), ao(980) [H, [3, i, i, 
H, @ and A(1405) 0, i, Ahhough it is still diffi- 
cult to identify an exotic state as a hadronic molecule, 
the exploration for possible molecules in more systems 
is an interesting topic. Such a study may help us to 
understand the strong interactions. There were dynam- 
ical studies whether the possible molecules exist in the 
light quark systems. In Refs. [l^, El, [H, various 
meson-baryon systems were investigated, while in Refs. 
[S d, m, [O, El, El, iOl , the nN and NN systems were 
studied. 

For heavy quark systems, the formation of molecules is 
easier due to the relatively small kinetic term. The rele- 
vant study can be traced back to thirty years ago [2ll[23|. 
Ten years later, Tornqvist studied possible deuteron-like 
meson-meson states bound by pions in Refs. [H, [23| 
which were called deusons [l^. In Ref. [1^, Ericson 
and Karl investigated the critical mass for molecule for- 
mation. In recent years, the renaissance of hadron spec- 
troscopy, es p eciall y the observation of exotic heavy quark 
mesons [13,111, [H^, , triggered extensive discussions in 
the molecular picture. 

The charmed meson Dsj{'2.i\l) [Sll, [s^l whose mass is 
much smaller than the quark model prediction was once 
proposed as a DK molecule [s^. Similarly, Dsj(2460) 
[33 | was suggested as a D*K state. However, their cs 
nature is strongly favored after considering the significant 
contributions from the DK continuum [3J|. 

The discovery of X(3872) [H, [H, [13, [11 ignited physi- 
cists' great interests. It is almost on the threshold of 
D^D*^ and very close to the thresholds of pJ loJ jij] 
and D~^D*~. The most popular interpretation for this 
intriguingstate is a hadronic molecule dominated with 
D"D*° [3i, [H, [Hi, [13, [13. However, this picture was 



questioned in Ref. [4J|. Very recently, the BaBar col- 
laboration measured a relatively large branching fraction 
for X(3872) '0(2S')7, which indicat_es that X(3872) is 
possibly a mixing state of cc and D^D*^ [isl ]. 

For the interpretations of the exotic Y(4260) [H, [13, 
lisj l in the molecular picture, Liu et al. suggested it is a 
XciP state [i^ while Yuan et al. proposed it is a Xci^ 
state [13] • There are also other molecular proposals such 
as a Ac pair [EH and a DqD* or DiD bound state 52, 5^. 
In fact, the most popular opinion is that Y(4260) is a hy- 
although this interpretation is also 
58j |. We still require detailed investiga- 
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inconclusive 

tions to answer whether these interpretations are correct 
or not. 

Recently, the Belle collaboration observed a charged 
charmonium-like state ^"''(4430) in the n'^Tp' invariant 
mass distribution [soj . This state is an excellent can- 
didate of heavy quark molecules. The dynamical cal- 
culation also indicates Z"*" (4430) may be interpreted as 
a DiD* {D[D*) molecule ^6^, [O, [gl, [H, [H . Not long 
ago, the Belle collaboration announced two more charged 
charmonium-like resonances ^"'"(4051) and ^+(4248) in 
the TT+Xci mass distribution [65,], which gives us the hope 
that heavy quark molecules do exist. Unfortunately, the 
BaBar data do not support the existence of ^+(4430) 
[11]. Cross-checks for the other two charged resonances 
are also desired. 

Therefore, none of the heavy quark molecules has been 
established yet. With the development of experimental 
measurements, more and more exotic states in the heavy 
quark region will be found. It is worthwhile to study in 
which systems heavy molecules can exist. Motivated by 
the observation of new exotic states and the possibility 
of forming heavy quark molecules, Wong explored the 
combinations of heavy mesons and heavy antimesons in a 
quark-based model and found many molecular states [isj . 
Voloshin and Dubynskiy suggested the possible resonance 
at the D*D* threshold in Refs. [13, Zhang et al. 
studied possible S-wave bound states of two pseudoscalar 
mesons using the vector- meson-exchange potential f6£j]. 
In Ref. [t^], a DgD* molecule was proposed. 

In a previous work [7l| , we studied the S-wave DD/BB 
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D*D*/B*B* and D*D/B*B systems in a meson ex- 
change model at hadron level, where we considered 
scalar, pseudoscalar and vector mesons exchanges. In 
this paper, we will explore similar systems in a chiral 
SU(3) quark model (xQM) 72] and calculate the binding 
energies by solving the resonating group method (RGM) 
equation [7^. All the mesons below 1.1 GeV will be con- 
sidered. The study can be used to test different model 
approaches. 

The chiral quark model is a useful tool in connecting 
the QCD theory and the experimental observables. It 
has been proved successful in studying the baryon-baryon 
interactions and the meson-baryon interactions. For the 
mechanism of the short range quark-quark interaction, it 
is still controversial whether one-gluon exchange (OGE) 
or vector-meson exchange dominates. In Ref. [7^], Dai et 
al. extended the chiral SU(3) quark model to include the 
vector meson exchange part and named the model the 
extended chiral SU(3) quark model (ExQM) which was 
also successful in reproducing the energies of the baryon 
states, the binding energy of the deuteron and the NN 
scattering phase shifts. 

It is interesting to study whether this phenomenolog- 
ical approach is applicable to the heavy quark systems. 
We have applied this model to the D^D*^ system in Ref. 
[tsI and we will continue to perform similar studies to 
other systems. One may test this model by comparing 
the predictions with future measurements. 

The paper is organized as follows. After the introduc- 
tion, we present a brief discussion about the systems we 
will study in Section [Til In Section IIIIl we present the 
ingredients of the model. Then in Section IIVI we give 
the essential parameters for the calculation. We show 
numerical results for different systems in Section |Vl The 
last section is the discussion and summary. 



II. HEAVY QUARK MESON-ANTIMESON 
SYSTEMS 



The S-wave single heavy quark mesons are pseu- 
doscalar type [£), Ds, B, Bg] and vector type [D* , 
D*, B*, B*]. For simplicity, P {V) wiU represent the 
heavy quark pseudoscalar (vector) meson. We investi- 
gate whether the hadronic molecules can be found in the 
combinations of these mesons and their antiparticles in 
this paper. From the flavor SU(3) symmetry, the mul- 
tiplets are 3x3 = 8 + 1. One may consult Ref. [Tlj 
for the explicit flavor wave functions. The largely broken 
SU(3) symmetry must be taken into account for possible 
hadronic molecules. In the numerical evaluation, we will 
first consider the isospin symmetric case. Because isospin 
symmetry breaking (ISB) is probably important, we will 
also discuss the case of large isospin breaking. 

In the isospin symmetric case, we need consider only 
four possibilities:_(l). 1=1/2 (Qm)-(Qs), _(2). 1=1 (Quj- 
{Qd), (3). 1=0 (Qs)-(Qs) and (4). 1=0 {Ql)-{Ql), where 
Q is a charm or bottom quark and I represents an up 



or down quark. We call them 1=1/2, 1=1, l=0(s) and 
I=0(/) states respectively in the following parts. 

When studying the possible heavy molecule composed 
of a heavy meson and an antimeson, we take a simple pic- 
ture where only color-singlet mesons are involved. The 
OGE and the confinement interactions occur inside the 
mesons, while the meson exchange interaction occurs be- 
tween light quarks of different mesons. To make the de- 
scription accurate, we label the heavy quarks 1 and 3 
while the light quarks 2 and 4. The quarks 1 and 2 
are bound to the meson and the quarks 3 and 4 to the 
antimeson. We do not consider the flavor-singlet me- 
son exchange between heavy quarks or between a heavy 
quark and a light quark in present investigation. The 
consideration is as follows. In the chiral quark model, 
the constituent mass of the light quark is related to the 
spontaneous vacuum breaking while the breaking gives 
small effects to the masses of the heavy quarks, which 
indicates the coupling of the sigma meson and the heavy 
quarks is weak. 



III. HAMILTONIAN 

The details of the chiral SU(3) quark model can be 
found in Refs. [tI, [tJ]. Here we just present essential 
constituents for the calculation. The Hamiltonian for 
the meson-antimeson system has the form 



Tg + V 



OGE 



V 



conf 



M 



(1) 



where Ti is the kinetic term of the ith quark or antiquark 
and Tg is the kinetic energy operator of the center of 
mass motion. M is the exchanged meson between light 
quarks. 

The potential of the OGE part reads 



1 



1 



1 



1 
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(2) 



where Fq = ^ for quarks and F| = — ^ for antiquarks. 
ruq (itiq) is the light (heavy) quark mass. The linear 
conflncment potential is 



yconf ^ _4pc . p 



Q V"9Q 



^qh) 



There are similar expressions for V^q^^ and V^^^ . 

For a molecule formed with {Qq) and [Qq) mesons, the 
light meson exchange occurs only between q and q. From 
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Refs. [lalTJ, one gets 



(a = 0,l,2,... ,8) (3) 



C{gch,'mT,^,A)- 



C{gchv,mp^,A)^^Xi{mp^,A,nj) + 



(4) 



617121714 

fchv m2 + m4 fchv ^ 2 'Tl2m4 



9chv 



xX2{mp^,A,rij )[<t(i) ■ CT{j)\][\S)Xa{j% 



M 



(5) 
(6) 



Where Gm is the G-parity of the exchanged meson and 

9ch A^TO 



C{gch,m,A) 
Xi{m,A,r) 
X2{m,A,r) 
Y{x) 



47r A2 



A. 



Yimr) YiAr), 

m 

Y{mr) - (—^ Y(Ar), 
\m J 



(7) 
(8) 
(9) 
(10) 



Here we do not present the tensor term and the spin- 
orbital term in the potentials since we consider only S- 
wave interactions. We use the same cutoff A for various 
mesons. Its value is around the scale of chiral symmetry 
breaking (~1 GeV). 

By solving the RGM equation, one gets the energy of 
the system and the relative motion wave function. From 
the definition of the binding energy, Eq = Mq^ + Mq^ — 
Msystem, one judges whether a system would be bound 
or not. 



IV. THE PARAMETERS 

There are several parameters in the Hamiltonian and 
the wave functions: the OGE coupling constants gq and 
(7q, the confinement strengths a^Q, the zero-point ener- 
gies a^Q, the quark masses mg and m,, the harmonic- 
oscillator width parameter 6„, the quark-meson coupling 
constants gch, gchv and fchv, the cutoff A and the mixing 
angle for the / = mesons. The mass of the phenomeno- 
logical a meson is also treated as an adjustable param- 
eter. For other meson masses, we use the experimental 
values. 

The sigma meson does not have a definite mass. In the 
light quark systems, this mass parameter was adjusted 
to fit the mass of the baryons, the binding energy of the 
deuteron and the NN phase shifts. When extending the 
application of this model to the heavy quark systems, we 



XQM ExQM 
Set 1 Set 2 Set 3 



bu (fm) 
mu (MeV) 
rus (MeV) 
(MeV) 

Qchv 
fchv I Qchv 



0.5 

313 
470 
595 



0.45 
313 
470 
535 



0.45 
313 
470 
547 



2.351 1.972 

2/3 



TABLE I: Three sets of model parameters. Other meson 
masses are: = 984.7 MeV, = 980 MeV, = 138 



MeV, rrir, = 547.8 MeV, m^/ 



957.8 MeV, = 775.8 MeV, 



= 782.6 MeV and = 1020 MeV. 



use the values determined in the light quark systems. If 
the vector meson exchanges are not included, the mass is 
595 MeV, while m„ = 535 MeV and 547 MeV were used 
in the ExQM. 

For the up and strange quark masses, we use the values 
given in the previous work Q, [zl, [z3' = ^13 MeV 
and nis = 470 MeV. To investigate the heavy quark mass 
dependence, we take severa l ty pical values rric = 1430 
MeV [Zi, mc = 1870 MeV (73, rub = 4720 MeV which 
is close to the value in Ref. [tJ, and TOh=5259 MeV [tJ- 

The chiral coupling constant gch is related to gNNir 
through 



9lh 



Att 25 Att 



9 9nn^ ml 



(11) 



with glfj^^/{4:Tr) — 13.67 determined experimentally. 
From this relation, one gets gch — 2.621. In the ex- 
tended SU(3) chiral quark model, one also needs the vec- 
tor coupling constants. We adopt two sets of the values 
used in the previous work, (gc/iu, /chi;)=(2.351,0.0), and 
(1.972,1.315) [ill. The corresponding sigma mass is also 
presented in Table HI One notes each set of parameters 
can reproduce the masses of the ground state baryons, 
the binding energy of the deuteron and the NN and YN 
scattering observables. 

The values of gq, gq, a'^q and a^g can be derived from 
the masses of the ground state baryons and the heavy 
mesons. The binding energy for a system of two color- 
singlet mesons is irrelevant to the internal potentials of 
the color-singlet meson because of the cancellation [t^ . 
Therefore these four values will not give effects to the 
final results of Eg. We do not present them here. 

Isoscalar states with the same J^'~'' will mix. The 
mixing angle for pseudoscalar mesons r/i and ryg, 0^^ , 
is taken to be —23°. Because the mixing angle 9^ for 
scalar mesons is still unclear and controversial, we use 
three values in the numerical evaluation: 0.0, 35.264° 
and —18°. The second number corresponds to the ideal 
mixing while the last one is taken from Ref. [t^. We 
use the ideal mixing angle 9^ = 35.264° for the vector 
mesons. In the scalar and pseudoscalar meson exchange 
potentials, we have adopted Aq = I where I is the unit 

matrix. To investigate its effects, we also use Aq = y/§I 
to calculate the binding energies. 
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To consider the dependence of the binding energy on 
the cutoff, we use two values A = 1100 MeV and A — 
1500 MeV. 



V. NUMERICAL RESULTS 

When performing the numerical evaluations, we calcu- 
late the binding energies with all possible combinations 
of the parameters presented in the previous section. Only 
when all the results for a system indicate it is unbound, 
we conclude the system is unbound. On the contrary, 
we say a molecule is possible only when all the results 
indicate the system is bound. 



A. PP systems 

The quantum numbers for the neutral states are 
J^^ = 0++. The pseudoscalar mesons do not exchange 
in such systems since the coupling of three pseudoscalar 
mesons is forbidden. 

For 1=1/2 states, we investigate D°D+, B+B° and 
B^Df. Such systems are possibly bound by only scalar 
mesons a and e. After solving the RGM equation, we 
find these systems are unbound with various parameters 
presented in the previous section. 



Isospin 
I = 1 




System 
D°D+ 


xQM 

X 


ExQM 

X 






B+B" 


X 


< 1.9 






B+D+ 


X 


X 


7 = 0(s) 




D7Dt 


X 


< 10.4 






-fc>s i>s 


< 13.3 


2.5 ~ 43.7 






B°Dt 


< 3.0 


< 22.6 


/ = 0(0 D- 


D+ + D''D'' 


< 4.9 


13.7 - 52.9 






B° + B+B- 


10.1 ~ 26.8 


47.2 ~ 102.3 






D+ + B+D° 


0.4 ~ 12.8 


23.1 - 72.8 


ISB 




D^D^ 


X 


< 12.8 






B+B- 


0.1 ~ 10.3 


11.7 ~ 43.3 






B+D° 


< 2.0 


0.5 ~ 24.2 



TABLE 11: The binding energies for PP states, x means the 
system is unbound. 

For 1=1 systems, we calculate the binding energies of 
D^D+ , B+B^ and B+D+. Vector mesons p and uj are 
permitted, but p exchange interaction is repulsive while 
UJ is attractive. Their contributions are almost canceled. 
a and e provide attractive force while a' gives repulsive 
interaction. By exploring different cases of parameters, 
we get the binding energies for these systems. The final 
results are given in Table [TTl If all the numerical values 
indicate the system is unbound, we mark it with "x". If 
the system is unbound with some parameters and bound 
with other parameters, we give the upper limit of the 
binding energy. If all the results indicate the system is 
bound, we collect the binding energies in a range. From 



the table, one knows there are no bound states in D^D+ 
and B+D+. 

The hidden strange l=0(s) states we investigate in- 
clude D~D+, B°B° and B°Df. They may be bound 
mainly by the attractive cr, e and (f>. According to our 
model calculation, it is difficult to draw a definite con- 
clusion whether the bound states may form (see Table 

m- 

The 1=0(0 systems we study are -^{D- D+ + D°D°), 
-^{B°B° + B+B-) and -^{B°D+ + B+D°). Comparing 
with the 1=1 systems, a' and p exchange interactions are 
both attractive now. The amplitudes of the potentials 
are also larger. From the results in Table HIl we find the 
bound states containing bottom quarks exist, even if only 
scalar mesons can exchange. 

In the real world, the isospin symmetry is also bro- 
ken. The mass difference between and is around 
5 MeV and it will affect the conclusion whether hadronic 
molecules exist or not. In this study, we also cal- 
culate preliminarily the extreme cases Z)°Z?°, B+B~ 
and B+D'^. Such cases get the minimum contributions 
from a' and p. Our results indicate the hidden bottom 
molecule BB is still possible. Table |IT] shows relevant 
results. 



B. VV systems 

The quantum numbers are J^^ = 0+^, 1^ , or 2++ 
for the neutral states. The pseudoscalar mesons, scalar 
mesons and vector mesons can all be exchanged in such 
systems. In our model, the amplitudes for scalar meson 
exchange interactions are the same as the PP case. 

Similar to the former case, we first investigate the 
D*°D*+, B*+Bf and B*+D*+ systems with 1=1/2. 
Here the vector meson exchanges are forbidden. The 
contributions from rj and rj' cancel largely and the pseu- 
doscalar mesons give finally small contributions. The a 
and e have not enough attractive force to bind the heavy 
mesons and these systems are unbound for the angular 
momentum J=0, 1, and 2. 

We explore three 1=1 systems D_*°D*+, B*+B*° and 
Comparing with 1=1 PP case, the exchanges 
of pseudoscalar mesons tt, rj and 77' are permitted. The 
contributions from rj and rj' reduce that from tt. For J = 
and J = 1, the interaction due to pseudoscalar mesons 
is attractive and for J = 2, it is repulsive. From the 
resulting binding energies, we conclude that D*°D*+ and 
B*+D*+ are not bound while B*+B*° is not excluded. 
We present our results in Table [1111 

The hidden strange states (1=0) include D*^D'*+, 
B*°B*" and B*^D*+. The contributions from 77 and rj' 
exchange interactions have the same sign. For J = and 
J = 1, they are repulsive. For J = 2, they are attrac- 
tive. Our numerical results are also presented in Table 
Hm D*-D*+ is not bound in xQM. 

For I=0(Z) systems, we calculate the binding energies 
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oi ^{D*-D*+ + D*"D*"), -1={B*"B*° + B*+B*-) and 
-^{B*°D*+ + B*+D*°). The tt, t] and V exchange in- 
teractions have hke sign. For J — and J — 1, they are 
repulsive while they are attractive for J — 2. We find 
there are no binding solutions for these systems in xQM 
if J = while the formation of molecules is possible if 
J — 2. Table Hill shows our results. 

Similar to the PP isospin breaking case, we study 
whether D*°D*°, B*+B*- and B*+D*" may be bound. 
According to our calculation, bound states in xQ^ do 
not exist if J = and the hidden bottom molecule is still 
possible if J = 2. We also present the results for this 
extreme case in Table Hill 



C. PV ± VP systems 

The components PV and VP do not have definite 
C-parity while the neutral PV ± VP states do. For a 
state with given C-parity, two conventions for the rela- 
tive sign have been used in the literature. The plus sign 
for the C=-|- DD* system corresponding to the X(3872) 
was widely used while the minus sign was adopted in 
Refs. do, Recently, Stancu analyzed the charge 

conjugation in multiquark systems in detail (8l| and 
she also obtained a minus sign. In fact, the conven- 
tion of the relative sign depends on the phase between 
P and P as well as V and V under the charge con- 
jugation transformation. But the final result is irrele- 
vant with the convention. For example, for the C=-|- 
D"D*^ state, one gets X = -1={D°D*^ - D*°D") with 

the convention D°{D*°) cu and D°{D*°) = uc. The 
resulting matrix element {X\(T2 ■ (J4,\X) is +1. If the 
convention D^{D*^) = cu and _D°(_D*'^) = cu is used, 
one gets X=^{D^D*" + D*^D^) and the same element 

{X\a-2 ■ a4,\X) = +1. In the following calculation, we 
adopt the latter convention which is consistent with the 
PDG assignment. So the quantum numbers for the neu- 
tral states are J^'-^ = 1+^ corresponding to PV ± VP. 

From the flavor SU(3) symmetry, it is easy to get the 
wave functions of other systems in the same multiplet. 
One may use PV + cVP to denote these wave functions 
where c is equivalent to the C-parity of the neutral state. 

In this pseudoscalar-vector case, the numerical results 
may be found in the PP systems or the VV systems. We 
explain this fact with 1=1/2 states. 

We investigate -^{D^D*+±D*^D+) and -^{B+Bf± 

B*+B'^y By comparing the binding energies with 1=1/2 
VV case, one finds the results for the c = +l(c=— 1) 
states are the same as those for J = 2 (J = 1) D*'^D*^ 
or B*^B*'^. Therefore these systems are also unbound. 
It is unnecessary to consider ± B*'^Df) 

since the mass difference between B^D*'^ and B*^Df is 
around 100 MeV and their mixing should be very small. 
For the system or B*~^Df, the results are the 

same as B~^Df of the PP case. 



Similarly,_for the 1=1 case, the results for the c = +1 
(c = -1)_D0D*+ and B+B*° are the same as J = 2 
(J = 1) D*°D+ and B*+B°, respectively. The results 
for B'^D*^ or B*^D^ are the same as B^D^ case. One 
can also get the results for 1=0 cases and large ISB cases 
from ,7 = 2 (,/ = 1) VV or PP. The correspondence for 
the numerical results between PV ± VP and VV or PP 
may be found in Table ITVl 

Such a feature is not difficult to understand. The dif- 
ference between the VV case and the PV ± VP case 
comes from the spin-spin parts of the potentials. The 
matrix element for the PV ± VP is {cr ■ cr) = ±1 while 
that for the VV is {cr-cr) = —2, —1 and -1-1 corresponding 
to J = 0, J = 1 and J = 2, respectively. Therefore the 
results for the c = -|-l(c=— 1) PV case are similar to 
those for the J = 2 ( J = 1) VV case. If pseudoscalar me- 
son exchanges are forbidden, the results for the PV±VP 
systems are similar to those for the PP. 



VI. DISCUSSIONS AND CONCLUSION 

From the numerical results in the previous section, we 
know the binding energy is always larger in the extended 
chiral quark model than in the chiral quark model. This 
is partly because the vector mesons provide attractive 
and relatively important interactions. Another reason 
is that the sigma mass in ExQM is smaller than that 
in xQM. This makes the attraction from a stronger and 
thus the binding energy is larger even if the contributions 
from vector mesons can be canceled. 

In order to make a clearer picture for the possible 
hadronic molecules, we summarize our conclusions in Ta- 
bles |Vl |Vl| and |VlTl In those tables, "x" means that a 
bound state does not exist. means a bound state does 
not exist in xQM while it is possible or not excluded in 
ExQM. "?" means we cannot draw a conclusion even in 
xQM and the system needs further study. "/" means a 
bound state is possible. 

From the tables, we know that the 1=1/2, 1=1 charm- 
anticharm and the 1=1 bottom-charm hadronic molecules 
do not exist. Therefore our conclusion for the D*~Df ± 
D~ Dl'^ system is inconsistent with Ref. 70]. Our cal- 
culation also indicates that the resonance-like structure 
Z+(4051) in the tt+XcI invariant mass [6^ could not be 
an S-wave D*D* molecule. 

On the other hand, the isoscalar hidden bottom 
molecules BB, J=2 B* B* , and C=+ BB* are very hkely 
to form regardless of whether the isospin symmetry is 
largely violated or not. All these states should be rather 
stable since B is the lowest bottom meson and B* do not 
decay via strong interaction. The experimental search 
for these states may be used to test our model. 

There are so many systems we cannot draw a conclu- 
sion, most of which are 1=0 states. Whether the effects 
due to coupled channels, the annihilation and the possi- 
ble mixing between S-wave and D-wave interactions may 
help is an open question. More detailed studies are nec- 
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Isospin System 
7=1 

B*+D*+ 


XQM 

J = J = l J = 2 

XX X 

< 3.5 < 1.2 X 

XX X 


ExQM 

J=0 J=l J=2 

XXX 

< 7.9 < 4.7 X 

XXX 


i = o{s) Dro:^ 

7-1*0 r)*0 

BfD:+ 


XX X 

< 10.4 < 11.8 < 15.0 

< 1.8 < 2.3 < 3.8 


< 10.2 < 10.3 < 10.4 
3.3 ~ 43.7 2.9 ~ 43.7 2.1 ~ 43.7 

< 22.5 < 22.6 < 22.6 


I = 0(Z) D''-D''+ + D'^D*" 
B "^B + B ^B 
B*°D*+ + B*+D*° 


X X 3.6 - 22.4 
X < 5.8 34.5 ~ 59.2 
X <0.4 13.9 ~ 36.7 


< 25.6 3.5 ~ 38.7 17.7 - 67.9 
22.6 ~ 66.4 34.7 ~ 84.2 60.2 ~ 120.8 
5.6 ~ 41.4 13.6 ~ 56.7 33.4 ~ 89.5 


ISB 7)*"7)*" 

B'+B*- 
B'+D*" 


X X < 1.1 
X < 3.5 5.3 ~ 20.2 
X X < 7.5 


< 6.7 < 9.6 < 16.2 
5.6 ~ 32.7 8.6 ~ 37.9 14.9 ~ 48.7 

< 16.0 < 20.0 2.3 ~ 28.6 



TABLE III: The binding energies for VV states, x means the system is unbound. 



Isospin 



1=1 



7= 1 
7 = 



PV + cVP 



7)07)*+ + cD*''D+ 
B+Bt° + cB*+B° 
DOD*+ + cD*''D+ 
B+ B*° + cB''+ B° 
D-DI+ +cDl-D+ 

(D°D*° + cD'-°D°)+{D-D*+ + cD*-D+) 
{B+B*- + cB*+B-)+{B°B*° + cB*°B°) 
B*+Dt/B+D*+ 
B*+7)+/B+7)*+ 
BfDt/B°D''+ 
(B*°D+ + B*+D'-')/{B°D*+ + B+D*°) 



c = +l 



7)*"7)J+ (J=2) 
B'+Bf (J=2) 
D*°D*+ (,J=2) 
(,J=2) 
DrDl+ (J=2) 
BfBf (J=2) 
D*°D*° +D'-D*-^ 
B*°B*° + B*+B*- 



c = -1 



7)*"7)*+ (J=l) 
B*+Bf (J=l) 
D*°D*+ (J=l) 
B*+B*" (J=l) 
Dt-Df (J=l) 
BfBf (J=l) 



(J 


=2) 


D*°D*° 




=2) 






B+ 


Dt 




B+ 


D+ 




B°Dt 



1 = \ 

7= 1 
7 = 



B''D+ + B+D° 

d"d: (i=i) 
B+m (1=1) 

D"D+ (1=1) 
B+B" (1=1) 
B+D° (ISB) 



ISB 



) 



D* 
B* 

^D°/B+D*° 
B+B*- +cB*+B- 



B* 



B*+B*- (ISB) (J=2) 



(ISB) (J=l) 
B*+B*- (ISB) (J=l) 



TABLE IV: The correspondence for the numerical results between PV + cVP and VV or PP. Here c means the C-parity of 
the neutral state of the multiplet. 



Isospin (c, c) (6, b) (b, c) 





X 


X 


X 


7 = 1 


X 


* 


X 


7 = 0(s) 


* 


? 


? 


7 = 0(0 


? 


/ 


/ 


ISB 


* 


/ 


? 



TABLE V: Summary of possible bound states in PP systems. 



(c,c) {b,b) {b,c) 



Isospin 


J=0 J 


=1 


J=2 


J=0 J 


=1 


J=2 


J= 


J=l 


J=2 




X 


X 


X 


X 


X 


X 


X 


X 


X 


I = 1 


X 


X 


X 


? 


? 


X 


X 


X 


X 


I = 0{s) 


* 


* 


* 


? 


? 


? 


■? 


? 


? 


I = 0(0 


* 


* 


/ 


* 


? 


/ 




'} 


/ 


ISB 


* 


* 


? 




? 


/ 


* 


* 


? 



TABLE VI: Summary of possible bound states in VV systems. 



essary. 

In the extended chiral SU(3) quark model, the gluon, 
pseudoscalar, scalar and vector mesons bind together the 
light quarks to baryons. The strength of OGE interaction 
in the SU(3) chiral quark model is greatly reduced due to 
the existence of vector mesons. But it is still controversial 
whether OGE or vector meson exchange dominates the 
short range quark-quark interaction. If the later mech- 
anism is not suitable, one finds more systems without 



binding solutions. 

In summary, we have performed a systematic study for 
the bound state problem of S-wave heavy quark meson- 
antimcson systems in a chiral quark model. The ex- 
changed mesons below 1.1 GeV have all been taken into 
account. Since we considered just color-singlet meson- 
meson configuration and several approximations were 
used, our investigation is preliminary. Our crude cal- 
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C 


= + 


C 


= - 


Isospin 


(c, c) 


{b,b) 


(c, c) 


(6,6) 




X 


X 


X 


X 


I = 1 


X 


X 


X 


? 




* 


? 


* 


? 


I = 0(0 




/ 


* 


? 


ISB 




/ 


* 


? 



TABLE VII: Summary of possible bound states in PV' ± 
systems. 



between heavy quarks or between a heavy quark and a 
hght quark is not excluded. Since no mass factor in the 
potential may suppress the sigma meson contributions, 
the value of the coupling constant ggga is crucial in dis- 
cussing whether or not such interactions are important. 
Although the coupling is expected to be weak, a small 
value may have big effects, which is also an open prob- 
lem in the present approach. 



culation disfavors the existence of 1=1/2, 1=1 charm- 
anticharm and 1=1 charm-bottom hadronic molecules 
but favors the existence of 1=0 BB, B*B* (J=2) and 
BB* (C=-|-) bound states. Whether the consideration of 
other effects, such as the coupling with hidden-color con- 
figuration and the coupling with possible D-wave, sup- 
ports these conclusions or not will be further studied. In 
our model, the sigma meson exchange interaction plays 
an important role in the bound state problem of the light 
quark systems. When extending the model to the heavy 
quark sector, the possibility of the sigma meson exchange 
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